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ABSTRACT 


Description  and  use  of  a  Fortran  general  sparse  solver,  modified  to  operate  effi¬ 
ciently  on  a  vector  processor,  is  given,  CRAY-1  performance  of  the  solver  in  analysis  of 
2-D  grids  is  presented. 


I.  INTRODUCTION 


A.  General  Sparse  Solvers 

The  performance  of  sparse  equation  solvers  on  vector  processors  is  highly  depen¬ 
dent  on  the  machine  instruction  set  and  associated  timings.  The  availability  of  indirect 
addressing  instructions  (gather /scatter)  is  a  particularly  important  issue.  When  the 
processor  has  a  hierarchial  or  other  flexible  memory  organization,  such  as  the  CRAY-1, 

control  of  the  data  flow  by  assembly  language  coding  may  also  become  important. 

Algorithmically,  another  set  of  considerations  is  introduced  by  the  possibility  of 
achieving  high  execution  rates  via  equation  re-ordering  and  recognition  of  favorable 
matrix  structures.  For  example,  blocking  of  matrices  can  reduce  data  traffic  [7]; 
recognition  of  global  sparsity  patterns  can  also  yield  a  vector  solution  [9]. 

B.  Generic  Vectorized  Solver 

In  view  of  these  opportunities  for  specialized  algorithms,  the  value  of  a  generic 
sparsity  algorithm  in  a  common  high-level  language  for  vector  processors  is  in  some 
question.  For  example,  it  is  known  that  traditional  general  sparsity  codes  [2]  execute 
poorly  from  a  high  level  language  on  current  vector  processors,  at  rates  less  than  1/200 
of  maximum  processor  performance  [10]. 

It  is  the  viewpoint  of  this  report  that  a  generic  vectorized  solver  has  value  in 
establishing  a  baseline  performance  against  which  the  performance  of  more  special¬ 
ized  solvers  can  be  judged 

The  following  properties  are  proposed  for  such  a  solver 

(a)  Consistent  with  general  solvers  for  scalar  processors,  the  solver  should 
accept  column-  cr  row-ordered  symbolic  matrix  descriptions,  plus  permuta¬ 
tion  vectors  describing  row  and  column  re-orderings. 


(b)  The  inner  loop  of  the  numeric  solution  should  recognize  opportunities  to 


exploit  vector  hardware. 


(c)  The  solution  may  be  of  a  two-phase  symbolic /numeric  nature,  where  a 
symbolic  phase  preprocesses  the  matrix  structure  and  provides  descriptors 
to  a  solution  phase,  for  possible  repeated  numerical  solutions  with  the  same 
matrix  structure  [2], 

11.  The  LU  Map  Approach 

A.  Introduction 

The  algorithm  to  be  discussed  is  a  variation  of  the  "looped  index"  or  "LU  map” 
method  of  Chuang  [l]  and  Gustavson  [2]  to  vector  processors.  It  was  first  discussed  in 
detail  in  1977  in  [3]  and  its  application  to  electrical  circuit  analysis  given  in  [4].  The 
following  discussion  is  taken  from  [3], 

B.  Symbolic  Vectorization 
1.  Introduction 

Given  a  matrix  A  and  ~ht  hand  side  B,  it  is  proposed  to  perform  a  triangular 
factorization  in  the  form 

A-LU  (1) 

where  L  and  U  are  lower  and  upper  triangular  factors  with  elements  /y  and  Uy ,  respec¬ 
tively  Column  ordered  reduction  is  performed,  with  /«  =  1.  The  forward  and  backward 
substitution  has  the  form 

LY-B  (2) 

UX-Y  (3) 

where  X  is  the  solution  vector  (For  descriptive  clarity,  pivoting  is  assumed  down  the 
main  diagonal) 


4 


I 


f 


i 
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2.  Scalar  Model 

In  the  work  of  Gustavson,  the  purpose  of  the  symbolic  phase  was  to  determine  the 
fill  charactristics  of  A,  i.e.,  the  exact  structure  of  L  and  U.  This  is  a  costly  process  that 
need  be  performed  only  once  for  a  given  matrix  structure.  To  acquaint  the  reader  with 
this  approach,  an  example  using  Gustavson's  "scalar”  map  is  shown  in  Table  1.  Special 
note  should  be  taken  of 

(1)  the  fill  positions  detected  by  the  symbolic  phase  in  the  generation  of  the 
LU  map; 

(2)  the  use  of  map  indices  in  the  numeric  solution  to  extract  information 
from  the  numeric  arrays  A,  L,  and  U; 

(3)  the  use  of  an  expanded  current  column  (X  array),  requiring  zeroing, 
expansion,  and  contraction  in  the  loading  and  storing  process; 

(4)  the  opportunities  for  the  use  of  vector  operations  in  the  numeric  solution, 
as  evidenced  by  the  indexed  array  operations  marked  "vector”. 

The  following  two  sections  are  intended  to  give  insight  into  the  symbolic  map  gen¬ 
eration  by  discussion  of  a  proposed  vectorized  data  structure  and  symbolic  operations 
on  it  during  the  factorization  process. 

3.  Vectorized  List  Data  Structures 

Consider  a  column  of  a  sparse  matrix  having  the  non-zero  row  positions  shown  in 
Figure  1  (before  fill)  This  structure  would  be  described  in  a  conventional  ordered  list 
as 

31,32,  .36,39,42,43,.. .47  (4) 

Such  a  list  enumerating  all  row  positions  will  be  termed  scalar  storage  Clearly,  the  list 
can  be  shortened  by  identifying  sets  of  contiguous  positions  (vectors)  and  retaining 
only  the  first  and  last  row  numbers,  viz, 
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(column-ordered  numeric  values  of  A  matrix) 

3. 2. 4  ,  2, 1,2, 6, 3, 1,3, 1,2, 3, 1,5 

(JA(j)  points  to  beginning  of  jth  column  of  A  in  IA) 
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(column-ordered  _lis_t  of  row  numbers  of  A) 
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(a)  Example  up  to  factorization  of  fourth  column 
Table  1.  Example  of  use  of  UJ  map  in  factorization 


loro  expanded  current  column  (X  array) 

Load  current  column  with  fourth  column  of  A 


indices 
from  IA 


X(2)-l 
X  ( 4 )  =  3~ 
X(5)  =  l. 

A 


vector 


Factorize  fourth  column 


indices 
from  IL 
of  pre¬ 
vious 
columns 


X(3)=X(3)-X(2)*L(2)=0-(1)  (1/2)  =  - 1  / 2l 
X(4)=X(4)  -  X(2)  *L(3)=3-  (1)  (1/4)  =11/4; 


X(4)=X(4) -X(3)*L(4)=ll/4- (-1/2) (*l/10)=27/10j 
X(S)=X(5)-X(3)*L(5)=l-(-l/2) (3/5) =13/10  j 

I  *  , 

;  starting 


indices 


from  JL 


01  ( 4 ) =1/X (4 ) =1 0/ 27 

X(5)=X(5)*DI(4)=13/27 
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indices 
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U(2)=X(2)1 

U(3)=X(3)J 

L(6)=X(5) 


vector 


indices  from 
IL, IU  of 


current  column 


(b)  Steps  in  Factorization  of  Fourth  Column 


vector 


Table  1.  Example  of  use  of  LU  map  in  factorization 


31,36,39,39,42,4? 


(5) 


This  form  is  natural  to  looping  operations  for  a  scalar  processor,  where  pairs  of 
numbers  are  directly  usable  as  upper  and  lower  loop  indices.  Alternatively,  the  initial 
row  position  and  the  vector  length  could  be  stored  as 

31,6,39,1,42,6  (6) 

This  form  is  favored  by  vector  processors  with  hardware  that  counts  down  vector  arith¬ 
metic  operations  to  terminate  a  vector  operation. 

Another  choice,  preferred  when  a  significant  number  of  singleton  (scalar)  posi¬ 
tions  are  present,  represents  a  vector  of  length  one  with  a  minus  sign  prefixing  the  row 
number  as 

31,36,-39.42,4?  (?) 

This  latter  structure  has  been  adopted  in  this  report. 

4.  Vector  Fills 

The  multiply-subtract  inner  loop  associated  with  factorization  can  result  in  pro¬ 
duction  of  fills  that  must  be  detected  in  the  symbolic  phase.  In  Figure  1,  the  process  of 
multiplying  the  kth  column  of  L  (termed  a  preceeding  or  recalled  column)  by  uk  r  arid 
subtracting  from  the  rth  column  of  L  (termed  the  current  column)  is  depicted  The 
zero-valued  positions  3?, 38, 40, 41,  which  initially  separate  two  vectors  and  a  scalar,  are 
filled  by  the  dense  vector  (36,43)  in  the  kth  column. 

The  symbolic  phase  produces  the  LU  map  by  scanning  the  numbered  pairs 
representing  the  vector  structure  of  all  the  preceeding  columns  and  the  current 
column  to  determine  zero-valued  regions  of  the  latter  covered  by  at  least  one  of  the 
former  These  are  the  fill  positions. 


Before  Fill 


After  Fill 


th  . 
r  column 


Figure  1.  Example  of  vector  fill,  with  data  structure 
description  (scalar  indicated  by  -  sign) 


III.  SOFTWARE  DESCRIPTION 


A.  Symbolic  preprocessing:  Generation  of  compressed  LU  maps 


CALL  NEWFOR(N,IA,JA,IPC,IPR.IPRI,JL,JU>J.IL,IU1JVA.JVL,JVU,IX,IPOS,LENSC) 
N*  is  number  of  equations 

IA(J)’  on  entry,  IA(J)  contains  row  number  of  jth 
column-ordered  matrix  element; 
on  exit,  IA  contains  compressed  vector-scalar 
format  of  same  row  number  information. 

JA(J)’  points  to  first  elements  in  IA  of  jth  column; 

JA(J)  is  changed  by  NEWFOR.  as  IA  changes;  JA(N  +  l) 

’User  suopiied  msut  data  to  subroutine. 


points  to  one  beyond  last  element  of  LA. 

IPC*  is  column  permutation  vector 

IPR*  is  row  permutation  vector. 

IPRI  inverse  row  permutation  vector. 

JL(J)  points  to  first  element  in  1L  and 

JU(J)  1U  of  jth  column;  dimensioned  at  least  N+l. 

1L(J)  contains  compressed  vector-scalar  row 

1U(J)  map  of  L  and  U 

JVA(J)  points  to  first  elements  in  numeric 

JVl'(J)  arrays  A,  U,  and  L  of  jth  column; 

JVL(J)  dimensioned  at  least  N. 

IX  scratch  arrays  of  dimension  N. 

IPOS 

LENSC  is  the  maximum  number  (&  1)  of  contiguous  non-zeros 
in  a  column  that  are  processed  in  scalar  mode. 


B.  Numeric  solution 

CALL  VMNP  (N,JA,IA,JVA,A,IU,IL,JU,JL,DI,U,L,X,JVU,JVL,IPC,IPRI) 

(see  above  argument  list  for  NEWFOR) 

A*  array  of  numerical  values  of  column-ordered  matrix 
D1  array  of  re-ordered  diagonal  elements  of  U. 

U  arrays  of  numerical  values  of 

L  column-ordered  upper  (U)  and  lower  (L)  triangular 
matrices,  diagonal  not  included 

X  scratch  array  of  dimension  N. 


C.  Forward  and  back  substitution. 


CALI.  VMBP(N,IU,IL,JU1JL.JVU1JVL.DI.U,L.B,X,IPC.IPR) 


'User  supplied  input  data  to  subroutine 


(see  above  argument  lists  for  1NEWF0R.VMNP) 

B*  array  of  numeric  values  of  right  hand  side  of  entry, 
and  of  solution  on  exit. 

IV.  PERFORMANCE 

Three  finite  difference  grids  [5]  illustrated  in  Figure  2  were  solved  using  this  code. 
The  equations  were  ordered  by  alternate  diagonals  j_6j,  which  yields  triangular  L U  fac¬ 
tors  of  the  form 

Di  Uiz 

Lgi  Ugg 

where  Dx  is  a  diagonal  matrix  and  Uzz  is  profile  matrix.  Although  Uz 2  can  be  solved 
more  efficiently  [7][8],  this  code  has  the  advantage  of  being  in  Fortran  and  simpler  to 
use. 

The  performance  on  the  CRAY-1  of  the  matrix  factorizations  step  is  depicted  in 
Table  2.  It  should  be  noted  that  the  execution  rate  is  approximately  proportional  to 
the  average  vector  length  during  solution.  Of  course,  this  is  not  true  asymptotically, 
since  the  rate  has  a  limiting  value. 


Problem 

Equations 

Time(msec) 

MFLOPS 

L 

#1 

391 

20 

1.59 

4.2 

#2 

507 

59 

5.39 

12.2 

#3 

2323 

652 

11.0 

27.0 

Table  2.  CRAY-1  factorization  performance  summary  for  three  grids  of  Figure  2; 


L  is  average  vector  length. 
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Appendix  A 
Program  Listing 


A  Z 

—  10 


Lc 


rc  £ 


intDr-cotnO’-CNfO'rinti) 


r*coo)  o  -  wn 

CM  CM  CM  CM  CM  CM 


• 

1 

. 

o  z 

CM 

1 

z 

3  < 

in 

CO 

1 

o 

X 

CM 

t 

1 

ft 

LU 

< 

CO 

GO 

1 

H 

a 

Z  -j 

CM 

cc 

o: 

CM 

1 

a 

z 

to 

< 

W 

UJ 

UJ 

1 

1 

DC  f. 

o 

z 

•  to  CJ 

CO 

m 

cn 

Q 

0) 

t 

o  ac 

UJ 

X 

3  a 

• 

z 

z 

Z  UJ 

o 

1 

u>  a 

X  >- 

H 

UJ  - 

3 

3 

<  H 

1 

3  co 

H-  -J 

►f 

Q  (O  < 

in 

z 

z 

< 

UJ 

1 

UJ  UJ 

z 

a 

Z  3 

CM 

o  z 

h- 

I 

Z  Q 

UJ  o 

o 

<  2 

^  CO 

UJ 

UJ 

UJ  f 

< 

1 

> 

a 

Q 

in  cm 

X 

X 

to  t- 

o 

J 

Z  X 

-1  o 

CJ  Z 

CM  W 

f  t— 

H- 

to  CO 

1 

<  ft 

O  UJ 

< 

Z  f  x 

CO  CJ 

+  CP 

UJ  UJ 

— • 

-f 

oc  a 

to  * 

-«  00 

z 

CM  z 

Z  CM  o 

o 

O' 

CO  *- 

o 

>- 

z  a 

~  3 

CO  \- 

I— 

a.  u 

o  < 

o  > 

H 

3  Z  = 

X  to 

•  CM 

z  ~ 

cc  Z 

z  z 

a 

fi 

CD  <  - 

UJ 

ft  - 

LU  f"  • 

O  H* 

a 

<  ft 

oa 

CO 

DC 

a 

■ 

CJ  J  UJ 

UJ 

CJ  to 

UJ 

z 

CC 

•  •-  Ul 

^  in 

•  3  a 

CJ 

ft 

<  co 

a  uj 

> 

< 

CJ  > 

in  cm 

UJ  Z 

a: 

O  00 

cc 

a. 

r-  -j 

z 

CM  CO 

z~ 

a  f< 

O  < 

i— 

a  uj 

to 

•  0)  o 

o 

CO  CM 

O  CP  CO 

to 

<  to 

z  a 

u  o 

z  > 

a.  to 

ft 

CM  w 

ft  CM  z 

z 

to 

ft  to 

<  UJ 

< 

-J 

(O 

W  oc 

CO  CO  o 

o 

-J  UJ 

lO  Q 
LU  Ul 

i-  oc 

LU 

oc  a 
o  or 

u-  O 

I 

Z  2 

<  z 
oc  3 
cj  -j 

o  o 
oc  u 
a 

< 

a 

UJ  on 
>  LU 
ft  <S) 
CC  lO 
Q  LU 
Ct 

<  CL 

s 

co  o 

«  u 

lO  K 

—  < 
X  X 


u_  uj  O  < 

z  a  to  a 

K  Z  UJ 

<  tt  I  uj  Z 

IO>l/)UJ 

a  C5 

< 

CL  Q 
t/)  Ui 
fS 


3  X 
a.  lj  h 
CO  Z  M 
Z  2 

> 

•  to 

o  ~  z 
z  >  o 

<  -I  f 
UJ  h- 

Q.  >  3 

z  •-»  -» 

Z  K  O 

>  U  to 

UJ 

to  0.  UJ 
Z  to  _l 

<  uj  a. 
a  a  m 

CO  ~ 

O  -J 

a  x  3 

a  ~  z 
a 

a  t-  a 

z  <  o 

<  Z  u. 


<  CC 

cc  a 
uj  »— 
Z  o 

UJ  UJ 

O  > 

a 

o  z 

h-  f 

CJ 

uj  Z 

>  o 

\  ft 

a  f— 

<  < 
-J  h- 

<  z 

U  LU 

to  Z 
3 
CJ 

o 

a 


•  >  UJ 

Z  a  to 

<  <  a 
X  3  < 

<  z  a 
-j  <  to 

<  3 

o  -I 
-  < 

•  Z  a 

<  <  LU 

o  z 

•  »-t  UJ 

a  x  CJ 

o 

=  HQ 

•  Z  LU 
UJ  N 
DC  U-  f 

a  a  a 
H-  O 
CO  •  f- 
3  CJ 
CJ  uj 
Z  -  >  ■ 

ft  CP  I 
*  0)  CC  I 

u  ^  o 

<  Li¬ 
en  h-  i 

a  —j  i 

x  o  <  i 

k  a.  3  i 
ft  uj  Z  i 
3  CC  < 
Z  i 


Z  CN  *-» 
UJ  W  to 
Z  3  Z 

f*  3  UJ 

o  -z 
^  in  ft 


~  z 
a  3 
a.  to 


ft  t— t  Q. 

CO 

Z  O  CO 
UJ  Z  f 
Z  <  I 


10  CM 
CJ  CM  CO 
UJ  Z  CO  CN 
>  3  CM  w 

<  (O  w  3 

I  •  u  > 
a  a  3 

CO  Z  ft  - 

>3 

<  co  ^  in 
a  -  in  cm 
a  x  cm  co 

<  -  co  CM 

3  CM  f* 

O  •  w  -I 
Z  -1  a  > 
f  -  CL  3 

3  f.  . 

o  a  ~ 
-j  •  z  in 

-l  CO  O  CM 
O  •  f  (0 
U.  <  to  CM 
Z  w 

ID  —I  LU  < 

X  <  z  > 

h  UJ  H  O 

cc  a  • 


~Q  3 

X  o  ~  z  ~ 


ft 

QC 

a 

u 

0. 

ft 

« 

>  —> 

3  a 

> 

*  a 

3 

3  f* 

• 

>  • 

< 

3  a 

> 

-  a 

3 

X  ft 

<  3 

-J  X 

z  - 

•  • 

3  00 

z  J 

•  • 

*  ft 

ft  -1 

CJ  • 

a  • 

Z  3 

*  3 

3  3 

—I  * 

(O  • 

3  f» 

•  J 

•  a 

a  3 

3  - 

z  > 

3  -J 

3  f 

•  > 

to  DC 

-1  3 

■  CL 

*f  * 

CJ  ft 

*  3 

X  <  X 
3 

CO  •  to 


<  I 

a 
a 
in  < 

CM  w  ^ 


co  a 
CM  z 


*-t  O  ui 

oc  O  > 

<  o 
x  o 

a. 

CL  * 

f»  cc 

3  > 
f  3 

H-  CM  < 

3  k  r- 

f 

f*> 

• 

CL 

< 

-  0. 

»  • 

<  X 

CP 

a 

< 

cc 

z 

z 

CC 

< 

Z  w 

a  to 

z 

Z 

CL 

CL  - 

-  3  Z 

<  to 

Z  CO  J 

• 

► 

f- 

f.  o  ^ 

<  -  - 

UJ  ft  > 

<  Ul  If  -“*>  « 

f 

- 

II 

II 

-  CL  CJ 

>  3  ~ 

X  *  < 

-i  •  O  O 

ii 

ii 

z 

3 

3 

<  ft  to 

3  3  n 

K-  ^  0C 

-i  ^  O  O 

3 

3 

• 

* 

• 

3  •  Z 

•  *  3 

O  OC 

X  O  O  O 

• 

- 

< 

fs 

*  <  UJ 

<  -1  • 

z  o  < 

u.  u  O  CP  CP 

ft 

3 

3 

<3-1 

ft  f' 

—  o 

O  3  ©  ao 

3 

3 

- 

<w 

«■ 

f  •  * 

*  *  3 

CM  CJ 

S  ip  - - 

w 

< 

< 

<L 

•  <  Z 

<  3 

to  Z 

CO  ~  3  -I 

< 

< 

3 

3 

f< 

CO  ~  - 

3  ft  CO  rr 

w  M  Q.  H  □  w  M  O  UJ  3' 


3 

ft 

3 

CC 

< 

oc  a 

CM  f  Z 

i 

w 

m 

z 

O 

O 

UJ 

o 

UJ  ft 

X  X 

<f>. 

/f- 

ft 

•  -I  z 

-l 

M  Z 

— i 

N  3  Z  LU  UJ  CD 

co 

<f 

CO 

ip 

f 

-1 

O 

-J 

i 

a 

_i 

1  O  O  W  CO 

co 

CL 

CO 

f 

II 

3 

3 

o 

*-t 

o  z 

o 

Z  UJ  ft 

-  f 

» 

z 

- 

w 

3 

ii 

n 

la¬ 

CO 

u. 

o 

CO 

li- 

o  a  to 

22  in  + 

in 

w 

in 

H- 

■f*. 

fv 

z 

z 

Z 

z  z 

o  o  ~  z 

< 

w 

< 

0) 

3 

3 

id 

UJ 

UJ 

UJ 

UJ 

>  UJ 

22  0  ii 

a 

3  Q 

z 

CO 

>— > 

•—> 

X  Z 

X  u- 

Z  X 

u-  <  Z 

Z  Z  <  **“ 

< 

ll 

< 

a 

CC 

O 

I- 

b- 

o 

if 

i— 

O  Z  f 

o  O  UJ  CL 

UJ 

< 

LU 

o 

o 

a 

CL 

Q 

a 

3 

u  cj  a  z 

CL 

z 

QC 

U- 

o 

ft 

ft 

Z  ^  i  -‘CJ 

uj  r*  — ■  »-• 

CJ  -  -  o 

Z  *  a  >cm 

<  —  co  Z  tp 


Z  uj  <  uj 
i  it  h-  3  K 


f*  CM 
H-  ll 
<  U 

z  to 

CC  Z 
O  uj 


•  w  CO  * 

<  a  a  z 

^Odw 

z  u-  f*  a 
oc  *  *  z 
O  uj 
U-  Z 


:  > 


<  < 

o  u 


3  -J 

>  -J 
3  < 

(J 


Z  ^  CM 

w  P-.  f 

o_  r-  uj 

as  •  in  co 

Z  r-  ~  - 

>  w  |- 

uj  <  O 
JKSh 

-»  a  o 

<  a  o  o  Z 

U  3u  O  u 


uuuuuuuuuuuouu 


CO 

CO 


0)  CO  f- 

GO  t-  f 

uuuuou 


■^MC0'Tincpr^c00)O*“CNC0’rincDr-c00)O*-cNC0tTLncpr“C00^O— cMCOtTincpr-cooo  —  cNCOrrcncpNcooiO’-CMCon 
T.^.<P-T.^^^^t-T.0j(Nr^(N(MCNCN(MrNcNn(onnnnnnnn^'t’to^’TTjTM''Tininininin 


a 

v 


<PAGE  1>  /////  FILE: SPARF  /////  >>>>>  MAIN  PROGRAM  <<<<<  ■ P/iGE 


END 


-*  a  2  cn  D^in(i)r-cooiO'^cNnTriDa)f^cocno  —  <n 

K'  cn  oo  — 


*.  o 

•  < 

a 


rn 


o 

a 


< 

z 


— *■  if) 
<J  Ui 
CL  H- 

•  Q 
QC  Ui 
CL  Z 
•-•  O 


<  if) 

•-I  o 

*  CL 

<  I 

“3  > 

w  -J 

z  z 

Ui  O 

z  z 

<  < 

oc  a 

Ui  (/) 

Z  ui 
>-*  ►— 
< 
3  a 
Q  ui 

a  z 

CO  Ui 
3  a 
00 


c/> 

~  o 

—  oc 

w  Ui 

<  fS 
"3  I 

-z 

~  O  I 

-  z  ; 

<  U_  l 

H  O  I 

^  a  : 

■*-  Ui  ' 

^  co  i 
u  z 
a  d 
~  z 

Ui 

-  o 

w  < 

a  of 

a.  ui 

—  > 
< 

z 

O  if) 
if) 

Z  * 

Ui  O 

Z  QC 
—  < 
a 


o 

oc 

Ui  Z  iO 
CO  -  ~ 
Z  J  ■ 
3  O  O 

z  a  • 

<  O 

Ui  — *  *“ 
I  n  ll 

in  h- 

i/>  *-•  < 

~  O  Z 

<  QC 
Z  Ui  O 

QC  u. 


O 

< 3 


O 
in  h- 


z 

< 


■  <3 


~  O  — *  - 

-  z  o  z  +  z 

I  -  z  ~  z  <  * 

*  *-  *  *  z  z  - 

O  m  *-  •  W  ^  ui  ii  it  “3  “3 

0)  o:  3  II  O  L  *■  it  3  '*•  3  II  II 

0)  <  »  JC  Ui  2  +  ^  Z  **  ^  ~  Z 

0)  CO  -<<<n  +  S330C 

o  II  Z  II  ’'0«-3Q:2Zi-  3Ww—3 
it  z  ii  -i  ^  wwiiwZ'-'  a  u  h  o 
<zzo<oou.u.<<o<oaauiz 

ZZ<l-30QMM2MU3QHHaUi 


n 

—  n 


(3 


—  cNOrrinipr-cooiO^wcorriniof-coOTO  —  MCo^ 
OOOOOOOOO’-**-*-  —  -  CNCMCiOiCN 


<PAGE  3>  /////  FILE; SPARF  /////  >>>>>  SUBROUTINE  RANGEN  <<<<<  <PAGE 


2  K 

2  Z 

3  Ui  UJ 

J2I 

O  UJ  h  > 

o  -j  a  ^  < 

mao  a 
a  z  k  <  a 
o  a  <  0  ^ 

u.  _J  uj  Z 

0  *  >  hm  co 

<  Z  ~~ 
M«aio  I 

W  O  H  UJ  h 

Z  »-  X  3 

hm  <  (j  -J  Z 

uj  <  *-* 
co  a  >  a  > 
a  o  a  a 

o  <  a  x  a 

J -3  m2 
Q-oasan 
*-»  a  o  o  j— 

a  a  a  <  uj 

0  O  Z  CO 
co  Z  Z  0  *-* 

iuhiz  <a 

Q  K  -I  *-*  CJ 


. 

* 

o 

* 

UJ  (/> 

* 

< 

< 

K  3Q  mi/1 

-J 

CO  CO  30  3  0  3-1 

< 

# 

K 

* 

X  tt 

0 

# 

Z  ►— 

UJ 

» 

Z  <  O  Z  UJ 

N 

UJIU20Z30Z 

# 

0 

* 

H  < 

Z 

* 

O  3 

X 

O  K  0  a  Z  Q 

■ 

ODoaujoctuiNH 

< 

# 

< 

• 

a 

HI 

* 

HM  Z 

»- 

•— »  t/>  Z 

_  oi  _  a  _  _ 

> 

« 

a 

* 

1/)  CO 

# 

f~  a 

w 

*-  OW  3<rt 

Z  Z  0  oo  0  co  z  Z 

3 

* 

* 

z 

o 

* 

— j 

<  UJ 

z 

< 

mI/)Zn  J> 

OOi^ZhiaZnZZ 

* 

3 

* 

a  a 

> 

# 

3 

*-  Q. 

* 

HM 

> 

1/»H  H  O  < 

3 

H  H  N  H  W  H  33 

< 

# 

-J 

* 

o  < 

HM 

# 

• 

3 

3 

O  J-  <  0  QC 

hm 

J—  H-  K  3  — J  -J  — l 

M 

# 

* 

u_  _J 

a 

* 

3 

S  J- 

CO 

• 

ai/lQ!w  Q1 

• 

mm  wa  3(/ia  30  o 

. 

# 

UJ 

« 

OC  3 

# 

3 

a  o 

HM 

z 

CO  <  h-  < 

co  co  m  <  cp  <  00 

< 

* 

co 

* 

UJ  0 

>• 

* 

- 

UJ  > 

01 

z 

Z<  I  U.  Ul 

oo<»-u.<ku- 

* 

a 

* 

UJ  CO 

OC 

* 

. 

h» 

. 

o 

HH 

-1  - 

•  *-  LL 

3 

3  3  -  “K  -  -  f- 

a  a 

* 

(O 

* 

Z  * 

K 

* 

t-M 

O  3 

0 

» 

O  UJ 

uj  z  a 

• 

OOujujZujujZ 

HM  HM 

♦ 

* 

HM  z 

HM 

# 

01 

>  O 

0  > 

>  UJ  O  3 

aa>>a>>a 

a  a 

* 

i 

* 

K  O 

00 

* 

a 

►m  a 

w 

U. 

M" 

HM 

Z  K  O 

■»“ 

HM  HM  Z  HM  HM  Z 

♦ 

* 

3  Z 

01 

» 

M 

a 

0 

o 

LL  H- 

h  UJ  01 

•w 

U.U.J-1-UJhMUJ 

a  a 

* 

* 

o 

< 

* 

• 

UJ 

a 

< 

o  < 

HM  U  01  0 

3 

00<Mj<Mj 

O  O 

* 

a 

* 

01  -1 

* 

0 

Z  co 

M 

CL 

3 

0  co  uj  uj 

3 

0  CO  UJ  0  CO  UJ 

— V 

* 

* 

03  < 

* 

# 

CL 

S  01 

UJ 

V-  UJ 

O  H-  uj 

M  t-  W  O  UJ  O 

X  X 

HM 

* 

z 

* 

3  a 

* 

* 

HM 

3  u i 

z 

CP 

z 

OC  Z 

a  »-  z  > 

z 

aaza»-zat- 

a  a 

01 

* 

# 

CO  UJ 

* 

* 

-J  > 

o 

z 

o 

< 

X  HM  O 

o 

<  <  X  X 

O  O 

a. 

* 

z 

* 

z 

» 

♦ 

01 

o  z 

HM 

3 

HM 

»-  a 

u.  uj  O  00 

HM 

KHU.LUJLlU.UJZZ 

HI 

# 

* 

CO  UJ 

« 

* 

UJ  -J 

0  HM 

CO 

z 

CO 

CO  *m 

mm  z  a  < 

CO 

C fl|/)MWZMMZ 

HM  HM 

- 

# 

> 

* 

HM  0 

* 

0 

Z 

z 

z 

*****♦##♦#** 

uuuuuuuouuuu 


w  < 

<  <  > 

“3  —  3  < 


u  u  u  u  u  u  u  * 


VD*J  “3  3*5  3-1  QX* 

* 

*■ 

uuuuuuuuouuu  0  0  0  0  000000 


ir)cor-coOTO“,“^r)'*tincDr^Go<7>O'^‘CMr5'tLr)CDt,^cocno,,~cNr)rrin0r,-QOO>O'^c>jr,3TTif)(Dr^cocj)O,"CNn'»3>!r>u)r'Coa^o^cN 

pjcNojoj<Nr)r)r)«nc*5r)oi)r)c,)n^^t,t^f^^r,T'?r'T,Tir)jrinininininminin0000U)0U)0U)0i>is-i^is*r'r-r^r^r^c^«oooo 


VECTOR  OF  UPPER  TRIANGULAR  NUMERICAL  VALUES'. 
VECTOR  OF  LOWER  TRIANGULAR  NUMERICAL  VALUES. 


(T)’Tini0Ma0)O**(Nnninu)N 

csNwwNcNOinnnoonnn 


*-  cm  co  in  ip 
xt  rr  O' 


r-  oo  O)  O  —  cn 

^  n  ^  in  in  in 


to  o 

“3 

z  o  a 

UJ 

a 

O  Z  n 

o 

w 

H- 

HH 

□  ID  H 

Z 

•-*  o 

3E 

CN 

< 

CO 

(J  < 

UJ 

3  U. 

3  O 

CO 

z 

o  z 

-I 

X 

O  Jh 

CN 

O  X 

CO 

N 

< 

o 

z 

3  »- 

—>  *t 

UJ 

t/y 

H  O  O 

CN 

K 

0 

< 

H 

3 

O  3 

n  uj  ^ 

o  cj  a 

o  o 

o 

"3 

O 

-J 

u_  O  3 

z  o 

UJ  -J 

<*>■  hh  i 

a 

o 

o  ♦- 

*- 

C3 

o 

O 

Z  h- 

X  3 

(J  ~ 

U. 

H 

UJ 

z 

o 

1 

< 

CJ 

~  z 

3  O 

K  3 

4  *-«  •*“ 

^Oi- 

O 

— ; 

^  M 

»- 

o 

"3 

HH  ►“  HH 

-i  a 

3 

3  O  ~  4 

h-  O 

C5 

3 

K 

i 

+ 

-J 

4 

ll 

K 

z 

o 

o 

z  • 

3  Z  -1  -1 

O  to  CD 

• 

o  < 

CL 

<  a 

O 

z  cd  a 

o  ^ 

*-•  3 

3  ^  -J  -J 

CD  Z  -1 

hh 

CD 

o 

< 

"3  « 

-1 

CO 

zoo 

o 

3 

-J 

1  w  w 

3  O  •-• 

o 

— J  ■ 

-l  u_ 

HH 

“3 

» 

"3 

< 

CL 

< 

in 

3  ~ 

z  ■ 

^  Z  3 

-J 

-j  a  -j  _» 

M  HH 

3  O 

d 

O 

wh 

—  »- 

i 

wh 

"3 

< 

M 

UJ 

-i  •  a 

^Ol- 

3  O 

■w 

-10  3  3 

•  h-  h- 

*  h- 

X 

cj 

*  -1 

*“» 

HH 

4 

4 

*3 

HH 

OhH 

3  a  a 

3 

3 

WH  II 

X 

4  ►— 

UJ  H  (|J 

w 

H  II 

H 

uj 

0. 

■w 

X 

11 

a 

< 

< 

1 

UJ 

a 

CJ  O  DC 

u.  • 

3 

3 

Oh  J 

3  CL  11  II 

-J  in 

-J  • 

-J 

3  ^ 

CO 

-l  *-• 

HI 

<  CL 

"3  Q. 

<  O’ 

< 

w 

-J 

K 

•  CJ 

3  3 

3 

3  CN 

H  H  J 

n 

o«ai- 

■  O  0. 

I— 

HW 

“3 

CN 

^  •  ' — ’ 

w 

*-H  < 

—  < 

HH 

“3  CO 

“3 

►H  H- 

z 

u.  a  i/i 

*-•  »“  3 

i 

<T 

—  tn 

3  a  o  o 

i  +  —  ^  -J  +  *  o  i 

'  X  "  CNXQT^-Q. 
uj  k  i  *-4  ; 


a  "3  +  a  +  ^ 

o  ii  ^  ^  ii  co  x  «  o»»  a  *  w 
'  ►-  a  -j  <  h-  <  «-i 


h  O  lii  d  Z  3  n  +  ii  O  "3  w  (j  I—  CD 

»—  OO'I  LUwii  J  O  "3  -I  H  IT  X  Si/m 

03  ^"31  V3i-  *3  UJ 


— 

11 

u_ 

< 

M 

— J 

-J 

U.  II 

-J  o 

o 

w 

M  LU  hh 

< 

<  u. 

II  <  o 

<  o 

<  w 

II  U. 

1- 

o 

H-» 

I  3  3  uj  u. 

3  3  0 

n 

3X0 

o 

X 

X 

HH 

• 

HH 

3 

3  O 

o 

X 

^  H  Z  i-  || 

H-4 

3  «-• 

^30 

3  O 

3  X 

X  ~ 

< 

'w' 

CJ  — 

3 

-1  HH 

3  3C5 

3 

3  O 

a 

6 

X 

z  z  < 

-1 

o 

< 

UJ 

< 

H 

2 

u- 

II 

:  3  < 

11 

* 

3 

* 

*- 

a. 

UJ  II 

UJ 

o 

Z 

-j  a  3 

X 

* 

CJ 

# 

HH 

3 

t- 

DC 

"" 

o 

« 

-1 

* 

< 3 

Z 

X 

Z 

-J 

CJ 

* 

< 

* 

z 

O 

UJ 

o 

WH 

UJ 

tn 

* 

(J 

* 

t— 

Z 

X 

1- 

O 

> 

* 

* 

X 

a 

a 

< 

O' 

CL 

* 

* 

h~ 

o 

H 

o 

(J 

UJ 

O 

* 

* 

o 

o 

UJ 

o 

o 

CN 

co 

ru 

O 

CN 

(o  (jo 

■«H- 

CN 

3 

* 

* 

z 

_i 

o  o 

3 

CN 

-J 

*— 

CN 

CN 

CN 

CN  CO 

CO 

CO 

3 

CO 

* 

•t 

o 

CJ 

o 

U 

o 

CJ 

CJ 

c 

c 

O  CJ 

CJ 

CJ 

o 

CJ  u  u 

c 

c 

c 

c 

o  u 

u 

CJ 

CJ 

CJ 

•“CNCO'tinipr-ooajO-'  otnrrif)wr*co0)O-~r>/nrTm&r^coo)O-c*n*TMV)^cDO>o  —  cNCo**iniDf**c00>O''-cNC0o’iniDr>'® 
'TTr'**rTrn"*^j’*Tinininininininininin<i>(DiDiDc0iDifiu>i0iDr*Hr'HfHrHrH.rwrHr'HrHf>HoocococDcococoffl<jocoa>a>O0>o>o>o>O)O,> 

OJC^CNCNCNCNCNCNCNCNrNCNCNCNCNCNCNCNrNCNCNCNrNCNrNrNPJCNCNCNCNCNCNrNrNCNCNCNCNCNCNrNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCNCN 


o 

< 

o. 


Z  O  ■-  cn  n 

in  a)  cn  cn 


rr  in  ip  r- 

O  0)  0)  0) 


I  4-  “3  o  UJ 

l  •  “3  ~  3  ~ 

•  O  Q  Z  O 

'll  II  U.  Z  M  CN 

i  «  “5  m  u  k  r* 

3  Z  Z  CN 

_»  ♦  o  a 

“3  *030- 

i  *  k  o  a 

— <  *  uj  h-  z 

x  *  a  oo  uj 


Nfl0®O**CNn^j) 

Lnininipuioioipio 

wnnnnnnnn 


CD 

* 

<  X 

* 

-J 

2  -J  -I  -  X  -1  -X  ■  • 

<  < 

2  < 

2 

* 

*- 

z 

* 

> 

22oaKoa»-ujuj 

>  > 

-J 

to 

* 

z 

o 

* 

_  Z  O  Z  O  2  2 

<  Z 

* 

O  D 

M 

* 

• 

zz^i-tt^i-au-i 

-J  -1 

>  o 

y 

# 

—  z 

to 

* 

mi 

nhWUOWUD<4< 

<  < 

u  * 

#  * 

*-  <  * 

a  * 

* 

■»“ 

UJ  U.  UJ  u.  >  > 

u  u 

1-  K 

< 

* 

< 

UJ 

* 

to  to  u_  >  u  > 

z  < 

a] 

* 

> 

* 

2 

aao  It  O  5  X  x 

a  a 

to  uj  2 

* 

a  a. 

* 

X 

> 

oo  «o  <  O  »-•  »-« 

UJ  UJ 

U4  Z  O 

CC 

o 

* 

UJ  to 

cc 

* 

»-« 

“3 

*-  »-  x  ax  ccce.cc 

z  z 

2  uj  uj 

o 

z 

* 

z 

UJ 

* 

cc 

• 

aiua  uj  a  b- 

2  2 

■J  -J 

b- 

< 

« 

X  Z 

o 

* 

H* 

«HQOOQOtJ<< 

w 

z  z 

<  uj  a 

o 

cc 

* 

«  > 

cc 

* 

< 

*— 

aazu.zzu.zzz 

M 

>  o 

UI 

o 

a 

* 

cc 

o 

* 

Z 

O  O  •"*  MM  M 

a 

CC  CC  Z 

2  U_ 

> 

K 

cc 

* 

l-  z 

* 

-I 

50  tO  5Q  5  O  2  -J 

<  < 

UJ  < 

U 

< 

* 

<  o 

<5 

* 

U4 

M 

UJUJ50Z50Z 

<*•>. 

-1  -1  X  o  z  cc 

Z  UJ 

3 

* 

Z  a 

Z 

* 

X 

• 

qooq:ujoq:ujhh 

2  2  1- 

MOO 

o 

> 

a 

• 

u. 

M 

# 

K 

cc  a 

OO  O  V)  <3  h 

M 

o 

* 

UJ  •— 

b- 

* 

— 

ZZ  C3  to  ciozz 

X 

z  z  z 

<  u 

1— 

z 

u_ 

* 

X 

o 

* 

z 

OOtflZNlOZnJZ 

- 

<  <  UJ 

Q  M  UJ 

< 

o 

* 

h-  X 

> 

# 

M 

2 

MM  22 

M  M  J  2  O  > 

K 

M 

UJ 

# 

— 

M 

• 

m 

Kf-  H  2  K-IU-J 

cc  cc 

< 

2 

1- 

to 

* 

to  a 

a 

* 

to 

■ 

MMioa2toa2oo 

W 

t —  1 —  u_ 

X  UJ  Z 

Z 

< 

DC 

* 

UJ  h- 

* 

cc 

z 

to  tO  CO  <  CD  <  uu 

m 

o 

to  o 

•W 

a 

H 

< 

* 

>  < 

> 

* 

a 

z 

— 

00<»—  u.  <  l-  u. 

. 

a  a 

i—  a  *-• 

a 

UJ 

2 

a 

# 

-J  Z 

DC 

* 

2 

a.  a  h(/)o«i/ioi-i- 

uj  uj  a  x  UJ  H- 

a 

a 

Z 

<o 

* 

o 

< 

* 

• 

-J 

-J 

to  to 

■*- 

a.  5  o 

O  >  2 

M 

a 

* 

<o  < 

a 

* 

u 

o 

2 

55  •  -i-  •  •  h-  a  a 

w 

0.  O  h- 

~  Z  -l 

. 

►- 

UJ 

i 

♦ 

H- 

* 

a 

o 

• 

OOujujZujujZmh 

-J 

2  -J  U 

a  o 

o 

Q. 

* 

uj  O 

M 

* 

M 

aa>>uj>>uju_u 

. 

UJ 

to 

> 

CL 

# 

Z  uj 

CD 

* 

- 

u. 

T“ 

M  M  J?  M  M  S 

u  u  > 

u.  u 

h- 

♦ 

a 

a 

# 

-J 

O 

W 

ItlLhhUJKI—  UJItlt 

*— 

o  a 

o  o  a 

o 

a 

O 

CO 

* 

»“  o 

< 

* 

“3 

2 

>— 

X 

UJ 

CL 

> 

*  a  2 

*  Ui  -J  Z 

*  ( 3 


i  o  to  uj 
uj  o  xx 
2Q.KUJLU 
X  o  Q 
I  IL  LL  Ui  ZZ 


a  cc  u  a  a  z 
o  o  h-  o  o  a 
h-  k  <  t-  h-  3 
UUQIUUI- 

Ui  ui  U  UJ  UJ  UJ 

>  >  to  >  >  a 


^  2  ^ 

~  x  a  — ' 

-Da  m  ~  +  * 

w  Q.  f-  t-  w  »-» , 

a  ►-  2  J  -  - 

zd.'—  i-  z  -j  -j  : 

*  *-«  CD  *-•  *  2  2  ' 

•r-  (- 

II  II  (/>  II  II  II 

II  X  ^  CD  *-■  Q-  1 

2  a  2  2  a  o  t 

d  •—  i/)  o  h-  co  ; 

in  m  x  n  *-•  m  • 

a 

o  a  o 

o  <  D 

3 
a 
o 


OU<JOOOO<JU<JUU 


iDp'Como  —  cNn^inur'Coo^O  —  <Nn-rinur-coaio  —  fNDrrmur^cocnO^cNOTyinor^coo^O'rcNn^inuJt^cocno  —  cnd 
o<oipur'r-r'p-r,^r-r'r-r-r~cDflooDcccDcDaocoooco<nai<T5a)a^aia)a)OT<^0000000000^-“^'^  —  —  ^'-•-'^cncncncn 
nnnnnnDnnnnnnnnnnnnnrtnnnnnnnnnnnnn’T^'fi^^n'f^^n^'rmr'JiT’i'jV'T'i 


cm  m  in  0  r-  coo)  O  cm  n  tj-  m  0  co  o  o  cm  n  in  0  r*-  co  cn  o  *-  cm  on  ><t  in  0 

—  —  —  CM  CM  CM  CM  CM  CM  CM  CM  CM  CM  D  (V  (V  D  ODD 


a 


CM  1 

i_i 

z 

A 

00  i 

X  • 

o 

a 

A 

0 

z 

Ui 

A 

CO  » 

X  3 

z 

O 

CM  ! 

f-  0 

□ 

z 

1  1 

»-«  < 

CJ 

o 

a 

0)  1 

3  > 

0 

Ui 

O  ■ 

z 

3 

cc  a 

M 

Ui 

O 

0  1 

o  o 

X 

-1 

K  i 

*-  h- 

> 

K 

<  1 

in 

<  < 

a 

0 

a  i 

(0 

Ui 

a  a 

0 

z 

i 

o 

<  < 

z 

z 

►H 

-  + 

Tt 

z 

a.  a. 

Ui 

> 

10 

< 

Ui  Ui 

o  a 

< 

• 

X 

in  in 

Ui 

Ui  o 

in 

n 

0 

Z 

►—  H- 

10 

a 

z  ui 

o 

<  in 

a 

Ui 

<  X 

a: 

Ui 

• 

f- 

X  h- 

f— 

<  - 

f- 

z 

H- 

in 

Q-  K 

u 

Z 

o 

NH 

o 

o 

ui  Z 

< 

U 

z 

a  t- 

z 

z 

in  ui 

< 

4—4 

> 

Ui 

< 

z 

Of 

X 

00 

*-  -i 

o 

Z 

y- 

Ui  0 

Ui 

Ui 

1— 

H- 

<  < 

< 

< 

a  uj 

Q- 

X 

o 

N-n 

Z 

ui  3 

Ui 

CJ 

<  in 

u 

*— 

z 

X 

H 

Ui 

oc  o 

3 

Ui 

3 

M 

Z 

a 

O  ui 

Ui 

> 

>-  a 

u. 

o 

w 

Ui 

a 

in 

3 

< 

ui  ui  o  x 

O 

^4. 

a 

< 

a 

3 

(/)  oc 

< 

X 

x  »- 

0 

H- 

V 

CM 

a 

a 

>-  o 

O  «n 

*- 

h-  a 

M 

s 

in 

r> 

H 

"“4  H 

3 

o 

>«✓ 

o 

o 

a 

in  < 

in  in 

0  0  3 

s 

z  z 

•w- 

z 

>  V 

0 

n 

in 

in 

X  z 

< 

i- 

2  X  »- 

\ 

X  X 

Ui 

a 

w 

X 

u  < 

u. 

o  a 

o  z 

ui  in  o  Z 

\ 

1-  h- 

O 

Z 

i-  — 

o 

Ui 

o 

o 

*-  X 

H  < 

hujOI 

0 

*-• 

tt  CC 

t— « 

uj  O  ui  \ 

UI 

O 

z 

o 

o 

3 

< 

0 

0  0 

1- 

>-<4 

Z 

QC 

z 

o  o 

• 

in 

r- 

a 

0 

0 

0 

in 

in 

0  O 

0.  *-» 

0 

o 

< 

0  0 

z 

X 

z 

>  in 

K 

~  < 

4^- 

in  i- 

<  3 

X 

0 

a 

< 

0  0 

CM 

in  o 

o 

4-H 

CM 

4—4 

•*“4 

>- 

4^ 

O'  0 

z 

Ou 

CL 

o  0 

NH 

in 

<  < 

cn 

1—4 

i— 

4-4 

+ 

\ 

X 

CL 

0 

z 

X 

a  x 

X 

a. 

0 

K  0 

0 

0  b- 

0 

X  0  0C  0 

1 

* 

X 

h- 

X  o 

0 

CM 

0 

a 

UI 

NH 

X 

0 

a. 

0 

NH 

X 

a 

0 

4-H 

X 

z 

z 

m  2 

in 

h- 

0  3  X 

OC 

X 

0 

a 

0  O 

X 

0 

0 

< 

—4 

X  \ 

z 

Ui 

in 

w 

a.  Z 

0 

in 

-s— ' 

a 

0 

in 

<w 

Q. 

0  in 

o 

4-4  0 

NH 

O 

—4 

>  Z 

NH 

0 

0 

NH 

0 

O 

Z  ui 

0 

OC 

w 

4-H 

UJ 

in 

4—4 

< 

0. 

UJ 

0 

in 

NH 

< 

a 

0 

in 

NH 

< 

CL 

0 

4-H 

4“ 4 

0 

•  z 

CO  X 

o  ~ 

0 

NH 

NH 

in 

nh  0 

NH 

z 

—4 

< 

w 

•f 

0 

X 

z 

4-H 

0 

in  0 

4—4 

0 

in 

4-H 

0 

in 

0 

0 

in 

0  0 

b- 

X 

K 

0  z 

II 

0 

> 

0  CJ 

X 

X 

0. 

< 

0 

4H 

0. 

o 

0 

II 

in 

4-H 

in 

<w 

0 

II 

in 

NH 

■w 

0 

II 

in 

4-H 

a  x 

X 

0  0 

0 

0 

4-H 

in  z 

0 

a 

II 

a  ui 

▼- 

UJ 

4-H 

+ 

UJ 

UJ 

< 

0 

<— «. 

4-H 

w 

< 

0 

4-4 

< 

0 

NH 

W 

<  0 

NH 

h- 

a 

in 

NH 

y-  o 

0 

in 

X 

UI 

o  -J 

w 

in 

in 

It 

w 

X  in  in 

X 

< 

k 

X 

< 

X 

< 

O  o.  »-« 

X 

■ - - 

H* 

0 

*— 

^4 

4- 

4”4 

0 

Z  0 

0 

in  0 

X 

Z 

in  -i 

< 

z 

0 

UJ 

V 

•— 1 1 

Nh 

4-H 

II 

<**■ 

0. 

II 

in 

II 

•—* 

a 

II 

II 

*“» 

a 

II 

in 

X 

X 

< 

NH 

a 

1 

X  oc 

0 

o 

0 

0 

NH 

y-  0 

0 

X 

nh 

C 

o 

Z  0 

< 

w 

>w4 

CC 

V 

0 

Ui 

fl 

0 

A 

X  UI 

(1 

V 

0 

0 

""4 

« 

3  3 

NH 

a  X  x  0 

0 

X 

a 

Z  0 

1 

Z  0 

in 

X 

0 

z  cj 

a 

4-H 

UJ 

Cl 

< 

UI 

0 

1-4 

in 

0 

CL 

X  X 

4—4 

in  x 

a 

4-H 

in 

0 

CL 

o 

O 

w 

II 

II 

3 

X 

-  3 

X 

NH 

0 

0  in  I 

0  in  2 

0 

X 

3 

o 

UJ 

in  uj 

II 

Z 

-J 

Ui 

0 

II 

II 

3 

NH 

>"— • 

w 

NH 

«-H 

0 

0 

4-H 

>W 

4-H 

4-H 

0 

4^ 

— H 

NH 

NH 

0 

♦— 

0  >* 

< 

•*4 

4-H 

w  O 

4-H 

o 

X 

0  < 

X  Z  Z 

X 

O 

UJ 

0 

Z  Nl 

in 

Z  UJ 

4—4 

in 

4-4 

a 

X 

II 

o 

0 

< 

W 

in 

X 

< 

w 

w 

in 

0 

< 

w 

—4 

in  o 

CC 

X 

0  z 

<  0 

< 

»— 

NH 

in  0 

4—4 

0 

w 

X 

NH 

a 

I-  CJ 

UJ 

UJ 

NH 

0 

-J 

■w 

4—4 

w 

UJ 

a 

a 

0 

|| 

4-H 

<w 

< 

< 

4—4 

X 

II 

NH 

— • 

< 

< 

NH 

4-4 

W 

< 

< 

4-H 

0  0  K 

II 

X 

X 

3 

II 

w 

II 

w 

W 

w 

0 

3 

0 

CO 

~  Z 

\— 

Z 

-J 

UJ 

II 

ii 

UJ 

< 

in 

UJ 

Ui 

0 

W 

I 

W 

W 

0  Z 

K 

NH 

4—4 

O  X 

CL 

0 

0 

X  in 

in  o 

NH 

II 

3 

0  UJ 

Z 

4—4 

< 

in 

-1  X 

0 

UJ 

4-H 

in 

in 

in 

X 

< 

u. 

in 

X 

< 

0 

< 

0 

Z  0 

H- 

4— 1 

o 

X 

u  3 

X 

0 

Z 

0 

H-  M 

0 

0 

X 

0 

in 

oc  cn 

«— * 

o 

NH 

Z 

►H 

nh 

O 

4—4 

in 

Ui 

4H 

4H 

4—4 

4—4 

NH 

NH 

NH 

NH 

4-H 

< 

4—4 

< 

< 

NH 

NH 

4-H 

NH 

4-H 

4-H 

O  X 

NH 

oc 

4-H 

NH 

4-H 

<  UJ 

0 

z 

O 

X 

0 

Ui 

0 

0 

X  K 

0 

0 

X 

Cfl  H- 

0 

0.  U 

NH 

o 

< 

UJ 

0 

)C 

X 

* 

0  OC 

> 

> 

y- 

3 

UJ 

Z 

0 

X 

Z  in 

X 

< 

X 

< 

< 

>:  < 

o 

o 

o 

0  Z 

o 

a 

1-4 

4-H 

O  0 

o 

4-H 

z 

NH 

o  z 

z 

0  CL  O 

z  o 

z 

*“ 

z 

NH  NH 

0 

-- 

CM 

n 

in 

in 

c 

c 

c 

0 

0 

(J 

0 

0  0 

0 

0 

0 

0 

0 

0  0 

0 

0 

0 

0  0 

0 

Nf-r-r-r-r'i'-r-cDoocoeieoeoeocoeocoaioiaioaiwoiOTaimoOOOOOOOOO--*  —  —  —  —  -  -’-'McvcNMWNtNiNcNn 


a  z  —ojcott  in  i£>  r-  co  <J)  O  wnTrmifii^cDaiO  ■^cMn-Tinor-cooO^cNn^Tin 
n  in  —  •'-•'-'-'-in  cNCNOJCNCNCNCNCNjcNDcncnnrom 


w 

in 

DC 

UJ 

— ‘  z 

3 

s 

<  3 

-1 

3 

Z  in 

< 

W' 

> 

+ 

< 

z  ^ 

+ 

•  »~ 

-1 

M 

U  W 

< 

-J 

2  O 

Z 

II 

o 

3  a. 

o 

U 

m  i-i 

a 

#-♦ 

•  • 

< 

M 

■w 

DC 

■w 

DC 

z  — 

o 

X 

z 

3  w 

1 

3 

in  a 

U- 

in 

in 

•  a 

u. 

» 

u  -* 

O 

*-s 

a  - 

z 

UJ 

3 

»-4  <•“** 

UJ 

o 

o 

<w 

-  ■*“ 

> 

M 

z 

< 

DC  w 

1- 

< 

+ 

0.  < 

H 

3 

z 

•-*  3 

< 

_J 

M 

»  • 

a 

o  ~  ~ 

z 

r>  w ■ 

< 

UJ 

in  3  — 

o 

■*-  u  — 

3  — 

Z 

w  + 

a 

o  +  z  + 

o  <  ~ 

K  H-l  3  *~» 

<  < 

a 

in  i  w 

UJ 

w  m  w 

UJ 

* 

o 

O  *  w  * 

• 

f- 

in  ^ 

a 

o  ~  ~ 

03  ^ 

3 

Z  3  0  3 

o 

~  3  X  3 

•  *** 

go 

O  W  U  w 

u. 

^  w  <  w 

<  w 

M 

-«  <  *4  < 

>-»  <  Z  < 

w  CD 

a 

f-  1  w  + 

o 

*"■  i  <  ♦ 

Z  ■ 

H- 

z 

<  ^  a  — * 

z 

w  -  #  *■*. 

DC  ~ 

in 

z 

3  *-  Z  -1 

< 

DC  -I  O  -J 

O  - 

i— • 

z 

O  ~  3  O 

CM 

a.  o  o  O 

u.  w  < 

o 

<  W  Z 

uj  —  (n  o  in  o 

in  z  ^  - 

-«  u  -  o 

<  z 

i 

Z  u.  * 

Z  i  ~  Z  "  — 

H- 

•  i 

UJ 

> 

*  Z  - 

UJ  •  ^  •—*  3  w 

o 

^  ■*— 

UJ  w 

z  z 

-j 

*-  <  »  o  o 

h-  —  •*-  —  3  in  _i  cc 

> 

II  <—»■  >-< 

3  Z  CO  4  CO 

«  O  II 

z 

ii  a  3  H  ii 

<  II  ^  +  H  W  II  O  II 

UJ 

>-H  l-«  —» ■  -f  || 

w  •  II  ^  II 

•  •  a  o)  o  i 

O  O  CD  ii  cn  3  3  n  3  ww  3  ►*  i 
I  ii  u.cDCD'-'O'*-''^-'— vZr>*<<  n  - 

^wh  cN'D’-au'Dam'DTrji 

«-  z  z  ~  z  z  ^  o  ii  »  o: 

I^D20<03DCflit0’'W0U: 

I  <  2Q  Q  to  1/1  QhhQm 


hh  ■:<  j^o^n 

•WW  n  o  J  •  -J 

<<«TiiuOnOO 

33n_J~CJ~U>- 

ii  li  o  W  —t  “3  *-1 

*“  CN  o  u  u.  w  ^  O 

^^O-*^Q3<C0C3 


UJ  UJ 

3  3 
Z  2  2 
~  a 
i-i-3 
Z  2  ^  O 
O  O  UJ  z 
UUOCUJ 


(DO'-ojDTriniDr-ooaiO'^cNrjrrinipr-cocDO^CNnniniDr-coaio^ojnTTiniD 

inipioiDioioioioioio(pr^f^r^e^Np-r^f-^i^a3oococooococooococoo)0)0)0<DW(D 

ifiinminininininininininininininininininminininininininininintninininininin 


PAGE  1 3>  lllll  FILE:  SPARF  Hill  »»>  SUBROUTINE  FORM  <<•:<<  <PAGE  13> 


r-i 


cm  d  nr  m  id  r-  ®0>O,-cMD'Tin<'>t*«.c00iO’,-cMr)^ini0r’‘C00iO''-'CMnnTin(Dr',*cc0>O*,"iMn,rin<0ri“C00iO‘,~<N 
’“’"’-’"'“’-^’-’■•-cjoicNWcNncNtNCNrtnnnnnnfonnnn'Tn’r^’TV'T’T^ininin 


ft 


i 

o 

■»“ 

i 

UI 

• 

- 

a 

• 

i 

/»■ 

UI 

CM 

i 

*- 

o 

00 

i 

W 

a 

1 

i 

- 

3 

o 

CO 

i 

-j  -j 

HH 

i 

CM 

i 

> 

• 

z 

1 

i 

a  o 

2 

01 

i 

z  * 

^, 

3 

c 

i 

<  < 

W 

-1 

i 

> 

-1 

O 

UJ 

►  “3 

o 

1— 

3  - 

- 

4 

-3 

UI 

2 

a 

4  — ■« 

■»- 

H 

a 

i 

- 

<w' 

< 

o 

4 

U.  -1 

3 

a 

u. 

O  hh 

-3 

UJ 

» 

z 

cr 

in  O 

UI 

< 

a  o 

*-*.  tx 

O 

-J 

<  • 

O 

< 

2  -i 

—1  w  1- 

O 

o 

“3 

3  a  o 

z 

in 

o  - 

-  a  lu 

< 

i 

UJ  b-4 

~  ~  > 

DC 

in  a 

X 

a 

in  a 

~~z 

UJ  M 

o 

a 

o 

a  * 

a^H 

K 

UJ 

a  oc 

a  in  k 

< 

> 

2  a 

a  < 

2 

o  « 

•  ah 

Z 

o  - 

^  —  o 

O 

o 

^  *  2 

UJ 

in  a 

a 

a 

UI 

UJ  >-« 

u  •*“  lu 

UJ 

QC 

K  - 

a  w  a 

c 

3 

<  < 

*-•  X 

a 

h- 

a  “3 

-  -  » 

o 

O 

UJ  « 

^  ^  .  o 

UI 

3 

z  4 

o  —  ^  a 

a 

a 

UJ  M 

in  w  ■*- 

i- 

a  ~ 

2  <  w  lu 

u. 

in 

oc 

36  "3  O  in 

O 

w  o 

-j  -  >  a 

3 

Z  U.  •  «  3  UJ 

-32-  ■  > 

•  Z 


3  Z 

O 

a  uj 

m  z 

3  ~ 

S)  H- 

3 

m  o 

—  DC 
X  m 
*-  x 
in 


<  -  *-  z  - 


a 

z 


-J  Q 
«  Z 
a.  < 


I  CM 


< 

•-I  x 


-  »-*  Dm 


•  Z  >  Z 

■0-3h 
~  *2 
u/)^a 

*  Z  m 
i  uj  ^ 

«  2  <  uj 
)  «  >  O 
Q  “3 
oO 


-  3 

II  w 

3  a 
a. 

o~z 

*— •  ft 

a  — 
o  Q.  0. 
o  ~  z 


nT 

o 


M  O  Z 
3  «  <-< 
^  2 
O  **3  a 
n  w  uj 
~  in  *- 
O  uj 

o  a  c 

Q  •- 


O  * 
n  ♦ 


-j  z  —  a  >-  —*  — •  d 

« ^  ^  -f  *-  o_  a 

U.  »■  "5  0.  a  •—  m  o  w  II  w  v 

O  ll  ■ — '  •— I  •— •  II  *•  II  l/)  x  H  X  X 

-DOww  •-.  *-* 

i—  OO  a  <  <  dd  ii  o  w  it 

(Si  O  O  «  «  *■«  vjow-'h  jz^i/ikOO 
*-•  II  II  >  >  11  II  II  11  II  —*  _J  H*  O  <  M  II 

_JDJO-JOa*-c'JDOOx<£oa2,-cM 
ZZZZO^«*-«*-iO^*-‘O^Oi-'^-l-J 


o 

o  >- 

I- 

o 

o  e> 

o 

~  *-  o 

-  £  • 

2  3  0 

3  •  UJ 

•  Z  • 

O  £  ~  ~  ' 


D 

O 


O 

e> 


CM 

¥ 


Z  »  D 

I  »-*  !-•  O 

3  2  Q. 

(I  W  D  *-• 
£  u.  O  u. 

~  a  ~ 


— *  + 

~  ~  o 


3  3  ¥  ~  +  i 

ii  n  w  ^  ^ 

n-  CM  U.  H  II 

1C  JC  M  ^  ^ 


in  — 

CM  -J 

O  O  -J  ' 


O  CM  O  <  < 

O  -i  0  2  2 

O  n-  *-  'I  CN  '  w 

-J  11  ^1£  OO 

‘t-  -  x  x  uj  uj 

J  <  3  3 

i  2  2  Z  Z 

_J  w  O  ^  II  11  >— I 

i  +  h  ip  ini-  x  x  h*  h* 

■  X  »  O  Ow<<ZZ 

II  CMO0.oau.22OO 
iSCjQnOnmhhUO 


a 

o 

in 

in 

z 

z 

UI 

UJ 

-< 

-/ 

3 

*- 

X 

i 

4 

3 

2 

z 

•• 

2 

4 

t-H 

3 

> 

> 

3 

-1 

z 

z 

3 

-1 

Z 

z 

3* 

-J 

> 

> 

3 

3 

in 

in 

o 

o 

a 

a 

3 

_j 

3 

3 

*  v- 

. 

3  ♦ 

X  -j  -I  ©  -j  »- 


o  ~ 

II  3 


3  — 
Z  w 
II  3 


in  -j 
o  -J 

a  < 

*-*  o 


o  < 

3  u 


o 

< 

a 


f- 

CO 

0) 

O^CMD^TiniOP^CO 

9 

0 

i 

CM 

D 

nr 

9 

s 

7 

8 

9 

o 

cm  n 

nr 

in 

10 

CO 

Gi  O 

•"  CM 

D 

nr 

in 

10 

r- 

a  G)  O 

CM 

D 

r? 

in  io 

CO  0)  o 

CM 

D 

nr 

0) 

0)  0) 

ooooooooo 

O  — 

—  -  —  CM 

CM 

CM  CM 

CM 

CM 

CM 

CM 

CM 

cm  n 

d  n 

D 

D 

D 

D 

D 

D  D  nr 

nr 

nr 

nr 

nr 

nr 

nr 

nr  nr  in  m  in  in 

in 

m 

in 

in 

lOlDlOlfilOlOlOlDlD 

9 

9 

9 

ID 

10 

10 

9 

9 

6 

6 

6 

6 

10 

O  (0 

10  10  10 

(0 

<o 

10  to 

10  10 

10 

10 

(O 

10 

10 

10  10  10 

10 

10 

10  10 

10 

10 

10 

10  10  10 

vO 

ID 

10 

10 

r 

L 


IU 


+  — -  + 


CL 

Ol 

-  . 

< 

t- 

w 

X  X 

“D 

< 

z 

■*- 

• 

<  < 

* 

*“ 

4* 

s 

• 

l 

■»- 

*“  1— ■  F-4 

< 

+ 

n 

-J 

o 

1-4 

II  ^  w 

< 

z 

V) 

II 

■*— 

o 

II 

-oo 

w 

z  o 

II 

UJ 

CO 

"D 

■f 

o 

CH 

X  Z  X 

3 

II 

CN 

3 

UJ 

“D 

“3 

o 

w 

<  1-4  < 

-1 

<■"» 

■— ^ 

•»- 

Z 

a 

0) 

w 

o 

00 

< 

~  z  z 

*— t 

4T- 

K 

Z 

♦ 

•-* 

a 

o  o  < 

< 

II 

11 

o 

Mri 

w  II  II 

0. 

< 

a 

K 

z 

o 

II 

“D 

X 

z 

II 

t/)  Z  x 

-i  Z  Z 

3 

■— * 

z 

o 

tl 

> 

li 

ii 

< 

1—4 

X 

o  *-  <c 

-i 

■w 

a 

l- 

o 

o 

< J 

< 

< 

o 

<N 

z 

z 

o 

< 

Q.  z  z 

< 

< 

o 

UJ  Z 

“D 

o 

z 

z 

a 

•-* 

*-* 

o 

CJ 

“D 

u- 

a 

UJ 

co  0) 
O  O 


<PAGE  15>  /////  FILE: SPARE  Hill  »>>>  SUBROUTINE  NEWFOR  <<<<<  <PAGE  15> 


